COMPLEXITY OF RANDOM SMOOTH FUNCTIONS ON COMPACT MANIFOLDS. I 



LIVIU I. NICOLAESCU 



ABSTRACT. We prove a universal equality relating the expected distribution of critical values of a 
random linear combination of eigenfunctions of the Laplacian on an arbitrary compact Riemann Tri- 
dimensional manifold to the expected distribution of eigenvalues of a (m + 1) x (m + 1) random 
symmetric Wigner matrix. We then prove a central limit theorem describing what happens to the 
expected distribution of critical values when the dimension of the manifold is very large. 
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Notations 

(i) For any set S we denote by |5| E Z>o U {oo} its cardinality. For any subset A of a set S we 
denote by I a its characteristic function 



I A : {0,1}, I A (s) 



i, s e A 
0, seS\A. 



(ii) For any point a; in a smooth manifold X we denote by 5 X the Dirac measure on X concen- 
trated at x. 

(iii) For any random variable £ we denote by E{£) and respectively var(^) its expectation and 
respectively its variance. 

(iv) For any finite dimensional real vector space V we denote by V v its dual, V y := Hom( V, M). 
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(v) For any Euclidean space V we denote by S(V) the space of symmetric linear operators 
V -> V. When V is the Euclidean space E m we set S m := S(E m ). 

(vi) For v > we denote by drf v the centered Gaussian measure on E with variance u, 

1 _ ;£ 2 

d7^(x) = e »«|da;|. 

Since lim^oT^ = we set 7o := <5o- F° r a rea l valued random variable X we write 
X G iV(0, v) if the probability distribution of X is 7„. 

(vii) If fx and z> are two finite measures on a common space X, then the notation fi oc v means 
that 

1 1 
-jl = v. 



1. Overview 

1.1. The setup. Suppose that (M, g) is a smooth, compact, connected Riemann manifold of dimen- 
sion m > 1. We denote by |dV^| the volume density on M induced by g. We assume that the metric 
is normalized so that 

vol ff (M) = 1. (*) 

For any u, v £ C°°(M) we denote by (u, v) g their L 2 inner product defined by the metric g. The 
L 2 -norm of a smooth function u is denoted by ||tt||. 

Let A g : C°°(M) — > C°°(M) denote the scalar Laplacian defined by the metric g. For L > we 

set 

U L = U L (M,g) := kcr(A-A s ), d(L) := dim J7 L . 

Ae[0,L 2 ] 

We equip U L with the Gaussian probability measure. 

T d(L) || u, || 2 

dj L (u) := (2tt) 2~e 2~|du|. 
Fix an orthonormal Hilbert basis (^ , fc)/ c >o of L 2 (M) consisting of eigenfunctions of A g , 

A^ fe = A fe ^fc, fc < ki => A fco < A fcl . 

Then 

J7 L = span{^ fc ; A fc < L 2 }. 
A random (with respect to d7 L ) function u £ t/ L can be viewed as a linear combination 

A fe <L2 

where u n are i.i.d. Gaussian random variables with mean and variance a 2 = 1. 

For any u £ C 1 (M) we denote by Cr(tt) C M the set of critical points of u and by D(u) the 
set of critical values 1 of u. If L is sufficiently large we can apply [15, Cor. 1.26] to conclude that a 
random function u £ U L is almost surely (a.s.) Morse, so that the random set Cr(tt) is a.s. finite. 

To a Morse function u on M we associate a Borel measure \x u on M and a Borel measure cr^ on 
E defined by the equalities 

\i u := ^2 S pi a u = u*(/j, u ) = ^2 n Cr(u)\S t . 

peCr(u) iSR 



^The set D{u) is sometime referred to as the discriminant set of u. 
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Observe that 

supp/Ou = Cr(it), supp cr u = D(u). 

When u G U L is not a Morse function we set 

fJ'u '■= \dV g \, cr u = 5o = the Dirac measure on R concentrated at the origin. 

Observe that for any Morse function u £ U L and any Borel subset B C R the number cr u (B) is 
equal to the number of critical values of u in B counted with multiplicity. We will refer to <r u as the 
variational complexity of u. 
We set 

/ t x rn_ / x 772 / ^ _ 771 

2 



Sm '- r(i + f ) ' dm 2r(2 + § ) ' ^ : " 4r(3 + § ) ' (U) 



(47r) _ T ^ (47r)"T ^ (47r) 

^'m • — ■ 

Let us observe that 

s m = h m (m + 2)(m + 4), d m = (m + 4)/i m . (1.2) 
The statistical significance of these numbers is described is Subsection 2.2. We only want to mention 
here that the Hormander-Weyl spectral estimates state that 

dim U L = s m L m + 0(L m ^) as L oo. (1.3) 

For L>0, the correspondence 

U L 3 u i-> fi u 

is a random measure on M called the empirical measure determined by the critical points of a random 
function. Its expectation is the measure p L on M denned by 



fdfi L = [ ( [ fdfiu) dj L (u), 
Ju L \Jm ) 



for any continuous function / : M — > R. Note that the number 



N L := [ dfi L = [ | Cr(u)\d-y L (u) 

Jm Ju l 



is the expected number of critical points of a random function in U L . 

In [ 1 6] we showed that there exists a universal constant C m that depends only on the dimension m 
such that 

N L ~ C m dim U L ~ C m s m L m as L — » oo, (1.4) 
and the normalized measures 

1 

N 1 

converges weakly to the metric volume measure \ dV g \ as L — > oo. This means that for L very large we 
expect the critical set of a random u G U L to be close to uniformly distributed on M. Additionally 
we showed that 

logC m ~ — to log m as m — > oo. 
Similarly, the random measure U L 3 u h-» cr u has an expectation 

cr L := EjjL(a u ) 
which is a probability measure on R defined by 



f{\)da L {\) = J (J f(\)drt(\)\ d~/ L (u) 



for any continuous and bounded function / : R — > R. Results of Adler-Taylor [ ] (see Subection 2.1) 
show that a L exists. 
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1.2. Statements of the main results. In this paper we investigate the statistical properties of the 
measure cr L as L — > oo and then as m — > oo. First let us point a small annoyance which we will turn 
to our advantage. 

Observe that if u : M — > R is a fixed Morse function and c is a constant, then 

Cr(c + u) = Cr(u), /x c+tt = /U u , 

but 

D{u + c) = c + D(u), cr u+c = 5 C * a u , 
where * denotes the convolution of two finite measures on R. More generally, if X is a scalar random 
variable with probability distribution ux, then the expected variational complexity of the random 
function X + u is the measure 

E(a x+U ) = v x * (Ta- 
in particular, if the distribution vx is a Gaussian, then the measure a u is uniquely determined by the 
measure E(ax+ U ) since the convolution with a Gaussian is an injective operation. It turns out that 
it is easier to understand the statistics of the variational complexity of the perturbation of a random 
u G U l by an independent Gaussian variable of cleverly chosen variance. 

Note that the lowest eigenf unction is the constant function 1. We consider random functions of 
the form 

A fc <L2 

where the Fourier coefficients Uk are i.i.d. standard Gaussians, and X u G iV(0, u) is a scalar random 
variable independent of the u^s. In applications u will depend on m and L. Equivalently, this means 
that we replace the Gaussian measure d-f L on U L with a Gaussian measure d'y^ of the form 

r 1 Kl 2 II"- 1 !! 2 

d it = — m , 6 2(1+0J) 2 \ du ^ 

(2vr) — VTTw 

where 

u = («o,*o) S) := « - ■"0*0- 
Since X w is independent of w we deduce that the expected variational complexity of X u + u is the 
measure <x^ on R given by 

at, = E(a x „ +U ) =7u*o- L . (1.5) 

Note that 



AT L = / doi(i) = / d<x L (t). 



The first goal of this paper is to investigate the behavior of the probability measures j^r cr^ as L — > oo 
for certain very special u/s. For reasons that will become clear during the proof we choose u of the 
form 

u = u{m, L, r) = Gi m ^ r L m (1-6) 
where r > and the quantity U3 m ^ r are uniquely determined by the equality 

d 2 

s m + Q m , r = r-^=:s? n . (1.7) 

Observe that as L — > oo we have u(m, L, r) — > oo so the random variable X u is more and more 
diffused. From (1.2) we deduce that 

™ + 4 - /r(m + 4) \ 

S m = r "^Sm. ^m,r = — 77" - 1 %■ (1-8) 

m + 2 \ m + 2) / 
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The inequality > s m shows that the parameter r must satisfy the m-dependent constraint 

m + 2 „ 
r > ——■ (Cm) 

771 + 4 

To formulate our main results we need to briefly recall some terminology from random matrix theory. 

For v 6 (0, oo) and N a positive integer we denote by GOE]y the space Sn of real, symmetric 
N x N matrices A equipped with a Gaussian measure such that the entries are independent, 
zero-mean, normal random variables with variances 

var(au) = 2v, var(aij) = v, VI < i < j < N. 

We denote by pN >v (\) the normalized correlation function of GOE^r. It is uniquely determined by 
the equality 

/ f(\)p N , v (\)d\ = ±-E GOE v N (trf(A)), 

for any continuous bounded function / : R — > R. The function pn,v(X) also has a probabilistic 
interpretation. For any Borel set B C R expected number of eigenvalues in B of a random A G 
GOE^ is equal to 

N [ PN>v (X)dX. 



Poo,v{X) :- 




B 

For any t > we denote by 3t t '• R — > R the rescaling map R B x h- > tx G R. If p is a Borel 
measure on R we denote by (Jl t )*P its pushforward via the rescaling map 3^. The celebrated Wigner 
semicircle theorem, [2, 13], states that as iV — > oo the rescaled probability measures 

(%^)JpN,v( X ) dX ) 
converge weakly to the semicircle measure given by the density 

1 

2ttv 

We can now state the main technical result of this paper. 

Theorem 1.1. Fix a positive real number satisfying r > 1. Let u = u)(m,L,r) be defined by the 
equalities (1.6) and (1.7). Then as L — > oo the rescaled measures 

converge weakly to a probability measure cr m r on R satisfying the equality 

r\ 2 

cr m , r oc 7 (r -i) * ( e 4> m+ i )r -i (A)dA ) , (1.9) 
where the symbol * denotes the convolution of two (finite) measures on R. 

Th above result has several interesting consequences. 
Corollary 1.2. As L — > oo the rescaled measures 

1 'a , W 



-/V V VdimC/i • 

converge weakly to a probability measure cr m on R uniquely determined by the convolution equation 

* cr m = tT m i, 

m + 2 V V"^ 7 * 



,\- 



where a mj i oc e * /j m+ i i i(A)dA. 
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Hence, the large L behavior of the average complexity <r L is independent of the background man- 
ifold M. 

Corollary 1.3. As m — > oo, the measures cr m converge weakly to the Gaussian measure 72. 

Let us briefly describe the principles hiding behind the above results. Theorem 1 . 1 follows from a 
Kac-Rice type formula of Adler-Taylor [1] aided by refined spectral estimates due to L. Hormander, 
[11], andX. Bin, [5]. Corollary 1.2 is a rather immediate consequence of Theorem 1.1 while Corollary 
1.3 follows from Corollary 1.2 via a refined version of Wigner's semicircle theorem. 

The basic facts coverning the Kac-Rice formula are presented in Subsection 2.1 while the proofs of 
the above results are presented in Subections 2.2, 2.3, 2.4. We have included two probabilistic appen- 
dices. In Appendix A we have collected a few basic facts about Gaussian measures used throughout 
the paper. In the more exotic Appendix B we discuss a family of symmetric random matrices and 
some of their properties needed in the main body of the paper. 

We want to comment on the similarities and differences between this paper and A. Auffinger's 
dissertation [3] which was a catalyst for the present research. 

Auffinger considers random fields on the round sphere S N with covariance kernel given by the 
function 

5i N>p : S N x S N -> R, K N , p (q, q') = \q • q'f, 
where • denotes the canonical inner product in D S N and p is a fixed real number p G 

{2} U [3, 00). Among many other things, Auffinger studies the behavior of the variational complexity 
of such a random function as N — > 00. 

The eigenvalues of the Laplacian on the round sphere S N of radius 1 are 

\ k (N) = k(k + N- 1), k > 0, 

with multiplicity 

, . 2k + N- 1 fk + N- 1\ k N - x 



If Lfc := yj Xk{N), then we can identify U Lk with the vector space consisting of the restrictions to 
S N C l^ 1 of the polynomials in (N + 1) variables and of degree < k. We have 

2k N 

dimU Lk ~ — - as k — > 00. 
N\ 

Using the classical addition theorem for harmonic polynomials, [14, §1.2, Thm.2], we deduce that the 
covariance kernel of the radom function denned by U Lk is 

k 

where denotes the "area" of S N and P n ^(t) denotes the Legendre polynomial of degree n and 
order £, i.e., 

P n/ (t) = {-lTR n (l) = 2-" r(n l + 2 / 1) (l - t 2 )"^ [j t j (1 - t 2 )" + ^. (1-10) 

In our paper we first let k go to infinity and then we let — > 00. In this case, our Corollary 1.2 is 
a statement concerning the distribution of critical values of the restriction to S N of a polynomial of 
large degree in (N + l)-variables. 

Let us point out that the limit k — > 00 leads to rather singular phenomena. The random function 
(field) defined by U Lk converges to a generalized random function a la Gelfand-Vilenkin, [10], whose 
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covariance kernel is the Dirac delta-distribution concentrated along the diagonal of S x S . The 
sample functions of this process are a.s. nondifferentiable so in in the limit the notion of critical point 
looses its meaning. 

2. Proofs 

2.1. A Kac-Rice type formula. As we have already mentioned, the key result behind Theorem 1.1 
is a Kac-Rice type result which we intend to discuss in some detail in this section. This result gives an 
explicit, yet quite complicated description of the measure <r^. More precisely, for any Borel subset 
BcR the Kac-Rice formula provides an integral representation of tT^(B) of the form 

*Z(B)= I f L ^B(p)\dV g ( P )\, 

J M 

for some integrable function /l,u>,b '■ M — > R. The core of the Kac-Rice formula is an explicit 
probabilistic description of the density f U L,B- 

Fix a point p G M. This determines three Gaussian random variables. 

(u L nL) 3 uu^uuip) gm, 

(U L ^t) B Uu) ^ duM G T;M, iBVJ) 
(U L ,^) 3« u 4 Hess P K,) G S(T„M), 

where Hess p (M w ) : T p M x T p M — > R is the Hessian of at p denned in terms of the Levi-Civita 
connection of g and then identified with a symmetric endomorphism of T p M using again the metric 
g. More concretely, if {x l )i<i< m are ^-normal coordinates at p, then 

m 

HesSpOOd^ = ''i., u 0J (p)d x i. 
1=1 

For L very large the map U L 3«4 du(p) G T*M is surjective which implies that the covariance 
form of the Gaussian random vector du LJ (p) is positive definite. We can identify it with a symmetric, 
positive definite linear operator 

S(duM) : T p M -> T P M. 

More concretely, if (x*)i<j< m are ^-normal coordinates at p, then we can identify S{ du^(p) ) with 
a m x m real symmetric matrix whose -entry is given by 

Sij(du u (p)) = E{d x% u 0J {p) ■ d xi u u {p)). 

Theorem 2.1. Fix a Borel subset BcR For any p G M define 

/l,w,b(p) := (det(27r < S , (ii w (p) ) )~ 5 | det Hess p (iiu>)| • Ib(u u (p) ) | du u (p) = V 

where E ( var | cons ) stands for the conditional expectation of the variable var g/ve« f/ze coti- 
straint cons. Then 

°t{B)= [ f L ^B(P)\dV g (p)\. (2.1) 
JM 

□ 
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This theorem is a special case of a general result of Adler-Taylor, [1 , Thm. 11.2.1]. The many 
technical assumptions in Adler-Taylor Theorem are trivially satisfied in this case. In [ 1 6] we proved 
this theorem in the case B = R and u = 0. The strategy used there can be modified to yield the more 
general Theorem 2.1. 

For the above theorem to be of any use we need to have some concrete information about the 
Gaussian random variables (RV^). All the relevant statistical invariants of these variables can be 
extracted from the covariance kernel of the random function m u . This is the function £^ :MxM-> 
K denned by the equality 

££(p, q) = E{u u {p) Ulli {q) ) = E((X + u(p) ) • (X + u(q) ) ) 

= u + J2 ^k(p)^k(q)=--^ + £ L (p,q). 

\ k <L 2 

The function £ L is the spectral function of the Laplacian, i.e., the Schwartz kernel of the orthogonal 
projection onto U L . Fortunately, a lot is known about the behavior of £. L as L — > oo, [5, 8, 11, 18]. 

2.2. Proof of Theorem 1.1. Fix L » 0. For any p £ M we have the centered Gaussian vector 

TO, w = o, 

{U L ,~/ L )3u^ (u(p),du{p),V 2 u(p)) eR®T*M®S(T p M). 

We fix normal coordinates (x*)i<i< m at p and we can identify the above Gaussian vector with the 
centered Gaussian vector 

(«(p), (d x iii(p)) l<i<mi 3 x i x j \ U {P) )l<i,j<m ) £ 8$M ® S. m . 

In [16, §3] we showed that the spectral estimates of Bin-H6rmander [5, 11] imply the following 
asymptotic estimates. 

Lemma 2.2. For any 1 < i,j,k,£ < mwe have the uniform in p asymptotic estimates as L — > oo 



E(u(p) 2 ) =s m L m (l + 0(L~ 1 )), (2.2a) 

E( d xl u(p)d x3 u(p) ) = d m L m+2 5 l3 ( 1 + 0(L~ l ) ) , (2.2b) 

E { d l^Mp)d 2 xkx Mp) ) = hnL^^Su + S ik d jt + SuSjk) ( 1 + 0(L- 1 ) ) , (2.2c) 

E{u(p)dl xJ u(P) ) = -dmL^dijil + OiL- 1 )), (2.2d) 

E(u(p)d xi u(p)) = 0(L m ), E{8 x Mp)0 2 xJx Mp)) = 0(L m+2 ), (2.2e) 

where the constants s m ,d m , h m are defined by (1.1). □ 

Now let u> = u>(m, L, r) be defined as in (1.6), (1.7). Using the notation (1.8) we deduce from the 
above that in the case of the random function we have the estimates 

E(u w (p) 2 ) =s% l L m (l + 0(L- 1 )), (2.3a) 

E( d xiUu> (p)d x juM ) = d rn L m+2 5 lJ ( 1 + OiL" 1 ) ) , (2.3b) 

E{d 2 xlxJ uMd 2 xkxl uM) = hmL^^dM + Sikd^ + Sud^l + OiL- 1 )), (2.3c) 

E{uMdl ixjUu} {p) ) = -d m L m+2 5 i3 {\ + 0{L~ 1 )), (2.3d) 

E{u„(jp)d xi u„(p)) = 0(L m ), E(d xi u w {p)dl xj u„(p)) = 0{L m+2 ). (2.3e) 
From the estimate (2.3b) we deduce that 

S(du UJ (p)) = d m L m+2 (l m + 0(L- 1 )), 
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so that 

i m 7n(m+2) . 

VI det S(u u (p))\ = {d m )^L—^ ( 1 + 0(L ) ) as L -> oo. (2.4) 
Consider the rescaled random vector 

(s L ,v L ,H L ) = (s L,LJ ' p , v L,u *'- p , H L,U,P ) := (L-fuUplL-^duM, L^V^M). 
Form the above we deduce the following uniform in p estimates as L — > oo. 

E{{s L f)=s^ + 0(L- 1 )), (2.5a) 
E(v^vf ) = d m S i:j ( 1 + 0(L~ l ) ) , (2.5b) 
E(H±HL ) = h m (S i:j S M + 5 ik 5 je + 5 ie 5 jk ) ( 1 + 0(L~ l ) ) , (2.5c) 
E{s L H^) =-d m 8 lJ {l + 0(L- 1 )), (2.5d) 
E{s L vf) =0(L~ 1 ), E{vfHf k )=0(L~ l ). (2.5e) 
The probability distribution of the variable s L is 

1 * 2 
dj s L(x) = =e 2a m(L) \dx\, 

v /27TS m (L) 

where 

Sm (L) = S - +0(L" 1 ). 

Fix a Borel set BcK. We have 

E(\detV 2 u u {p)\I B (u u {p) | d« u (p) = 0) =L^^,E;(|det J ?/ L |/ B ,(s L ) | v L = 0) 



L 2 / £ detff L s L = x, t; L = 0) ; -. dx . ( 2 -6) 



=:q L , p (L--r B) 

Using (2.4) and (2.6) we deduce from Theorem 2.1 that 

/ \ — 

*t{B) = L m U-L.J 2 Jj L , p (L-^B)p L ( P )\dV g (p)\, 

where pi : M — > E is a function that satisfies the uniform in p estimate 

p L (p) = 1 + OiL- 1 ) asL^oo. (2.7) 

Hence 

^ - (^y 2 £ <tt,(B)«<p)l«Bi<p)|. (2-8) 

To continue the computation we need to investigate the behavior of qi )P (B) as L — > oo. More 
concretely, we need to elucidate the nature of the Gaussian vector 

(H L | s L = x, v L = 0). 

We will achieve this via the regression formula (A.3). For simplicity we set 

Y L := (s L ,v L ) € l©l m . 

The components of Y are 

Y L = s L , Yt =vf, 1 < i < m. 
Using (2.5a), (2.5b) and (2.5e) we deduce that for any 1 < i, j < m we have 

E(Y L Yf) = + 0(L _1 ), E(Y^Yf) = d m 5 l3 + 0(L _1 ). 
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If S(Y L ) denotes the covariance operator of Y L , then we deduce that 

S ( Y Vi = + OiL- 1 ), S(Y% 1 = -Lfy + 0(L~ l ). (2.9) 

We now need to compute the covariance operator Cov(H L , Y L ). To do so we equip S m with the 
inner product 

(A,B) = tr(AB), A,Be§> m 
The space S m has a canonical orthonormal basis 

Eij, 1 < i < j < m, 

where 

E ={ Eij ' i = j 

13 i<j 

and denotes the symmetric matrix nonzero entries only at locations and and these 
entries are equal to 1. Thus a matrix A G § m can be written as 

A = a ij E ij = ^2 °*3 E Vi 

where 



V2dij, i < j. 



The covariance operator Cov(H L ,Y L ) is a linear map 

Cov{H L , Y L ) : M © M m — > S m 

given by 

Cov(H L ,Y L ) (f> a e a ) = £ E{H^)y a %j = £ E{H^)y a E ih 

\q=0 / i<j,& i<j,ct 

where eo, ei, . . . , e m denotes the canonical orthonormal basis in R © IR m . Using (2.5d) and (2.5e) 
we deduce that 

Cov(H L ,Y L ) (jry<* e <*) =-Vvd m l m + 0{L- 1 ). (2.10) 

\a=0 / 

We deduce that the transpose Cov(H £ ,Ye) v satisfies 

Cov{H L ,Y L Y \^h 3 Eij \ =-d m tr(A)e + O(L- 1 ). (2.11) 

The covariance operator of the random symmetric matrix 

Z L = z l,x ._ ( H L\ g L = x ^ y L _ ^ 

is then 

S(Z L ) = S{H L ) - Cov(H L , y L )5(y L )~ 1 Cov(H L ,Y L y. 

This means that 

E{zffi e ) = (E ij ,S(Z L )E ke ). 
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Using (2.9), (2.10) and (2.11) we deduce that 



CoviH^Y^SiY^CoviH^Y 1 ^ [ J^^Ai I = 4r tr(^)l m + 0^) 



E((zf-) 2 ) =h m + 0(L- 1 ), E{ 



z L -z L -" 

11 33 ' 



d 2 



h m -^ + 0(L~ 1 ), Vi<j, 

S !rr . 



and 



E{(zt) 2 )=3h m -^ + 0(L~ 1 ), Vi 
E(zM e ) = OiL- 1 ), Vi<j, fc<€, (i,j)^(M)- 



We can rewrite these equalities in the compact form 

d 2 



^(44) = ( K 



Note that with r defined as in (1.7) we have 



_ <4 (L2) r-1 



We set 



so that 



(r-1) 
2r 



E ( z tj z ke) = ^ K h m 5ij5 M + h m (5 ik 5j£ + + 0(L _1 

Using (A.4) we deduce that the expectation of Z L is 

rL\ _ n mi vi\o/vi\-l' x 



l m + 0(L- 1 ). (2.12) 



r(m + 4) 

We deduce that the Gaussian random matrix Z L " X converges uniformly in p > as L — > oo to the random 

X ■ 

r(m+4) " 



matrix ^4 / ^ l m , where A belongs to the Gaussian ensemble g^ hm > hm described in Appendix 



B. Thus 



lim qL, P (B) = qoo(B) := / E 2 Khm ,h m ( I det( A 



L— >oc 



B 



lm) |) 



r(m + 4) ' 1 ' yJ2-Ks 



dx 



(h m )~2 / £^ s 2«,i( detM 



lm) I)- 



(/if 



e 2 



2k, i ( I det ( A - a m yl m ) I ) —=dx, 
m v2vr 



where 



This proves that 



a, 



r(m + 4)\fh r 



lim ^ -1^(5) = (^m) 2 / Eg£*A 



L^oo (.3%) 2 



(1.2) 1 



detf A - -^=l r 



e 2 



=:/i m (B) 
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Using the last equality, the normalization (*) and the estimate (2.7) in (2.8) we conclude 



1 / h r , 



L^oc S m L m V (^^ m ) 2 



(1.2) / 2 \ 2 



(2.13) 



m + 4 , r^i + T )M«). 



Observe that the probabilty density of fi m is 

dUm _ 

, = £/ s 2k,1 

ay m 



det( A - -^=l r 



y 

e 2 



2tt 



(2.14) 



We now distingush two cases. 

Case 1. r > 1 From Lemma B.2 we deduce that 



E „2k.1 
dm 



det( A 7=l m 

\ IT 



m+3 _ /m + 3\ 1 

2— r 



(2.15) 



where 



Thus 



r 2 := 



k r — 1 



k — 1 r + 1 



±3/m + 3\ 1 (r 2 + l-2r)j; 2 



r/y 



2— r 



2 y 2vr\/K 



e 4r 



_ i (A-(r 2 +l)-^) 2 



/ m + 3 \ 1 



2 y 2^^ 

An elementary computation yields a pleasant surprise 

2 



/3m+1 , 1 (A) e -^( A -^ 2+1 )*) 2 -2(^TT ( iA. 



1 

'4^2 



A-(r 2 + l) 



Now set 



We deduce 



(A := y^A) 



V^J 2(r + l) = 
P = P(r) 



dUm _ ( m + 3 A L_ 

dy \ 2 y 2tt^ 



1 



-X 1 



2(r-l) a V 2(r-l) 



1 

-A - y* 



(r-1) 

p m+ i,i(A)e^ A2 e-K A -^) 2 dA 



m + 3 _ / m + 3\ 1 

2 2 r 



2 y 2ttVk 



f Vr~p m+ i,i(Vr~\)e- r ± x2 e-^ x -^ 2 d\ 

JR 



(B.6) m+3 _ / m + 3 \ 1 /" ... r \2 , , 

= 2— r — — / ^ A e-J* d 7 i 1/ " A dA. 

Using the last equality in (2.13) and then invoking the estimate (1.4) we obtain the case r > 1 of 
Theorem 1.1. 
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Case 2. r = 1. The proof of Theorem 1.1 in this case follows a similar pattern. Note first that in this 
case k = so invoking Lemma B.l we obtain the following counterpart of (2.15) 



det( A - yl 7 



m+4 ( m + 3\ 
2 2 r I — ^— ) e 4 p m+1> i(y). 



Using this in (2.14) we deduce immediately (1.9) in the case r = 1. This completes the proof of 
Theorem 1.1. 

□ 

2.3. Proof of Corollary 1.2. We use Theorem 1.1 with r = 1. Using (1.5), (1.6) and (1.8) we deduce 
that in this case 

= l 2s m L^ * a L . (2.16) 

m+2 



Using the equality 



m A- 1 

rm _ rm ~ f 

m + 4 



we deduce 



^(^),^ = (^),(^(^),^l- <2 - 17) 



Using (2.16) we deduce that 



Using the spectral estimates (1.3), the equality (2.17) and Theorem 1.1 we deduce 
-r fa i ) a L ) = lim — T 



lim 7 2 * ( -^f [H i_ ) <r L I = lim — -( 31 i ) af\ = (31 r—r) a m+1 1. 



We can now conclude by invoking Levy's continuity theorem [1 2, Thm. 15.23(ii)]. □ 

2.4. Proof of Corollary 1.3. By invoking Levy's continuity theorem and Corollary 1.2 we see that 
is suffices to show that the probability measures er m i converge weakly to the Gaussian measure 72. 
Set 

Rm{x) := y/mp m+ i tl (y/m x) = p m+1 j_(x), 

' m 

Roo{x) = WzI{\x\<2}V^ ~ X 2 . 



2ir 

Fix c G (0, 2). In [16, §4.2]. we proved that 



We deduce that 



lim sup \R m (x) — Roo(x)\ = 0, (2.18a) 

m ^°° \x\<c 

sup \R m (x) — Roo(x)\ = 0(1) as m — y 00. (2.18b) 

a 2 /4tt / A \ 1 a 2 

^i,i(A)e-r = ^-R m (-=) -=e" (2.19) 



and 

I fw f - I rn mx 2 4tt f - 

:= / /Wi,H A ) e ±d\ = \ — / R m (x)A — e * dx = \ — / R m (x)dj2_{x) 

Jr V m J R V 4vr Vmj R 

The estimates (2.18a), (2.18b) imply that 

I m ~ \/47ri2oo(0)m _ 2 as m — > 00. 
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To prove that the probability measures 

1 A 2 

-j-p m +i,i(X)e~~dX 
converges weakly to 72 it suffices to show that the finite measures 

1 A 2 

u m := m2p m+11 (\)e~~d\ 
converge weakly to the finite measure 

A 2 

Uoo := Roo(0)e~~dX. 
Let / : R — > R be a bounded continuous function. Using (2. 19) we deduce that 

/ f{X)du m (X) = f f(X)R m (m-h)e-^d\. 
Jr Jr 

We deduce that 

f(X)du m (X) - [ f(X)du 00 (X) = [ /(A) (R m (m--2x) - RooiO)) e~^dX 
Jr Jr v 7 

= / fWl{\\\<cy/j^ (-Rm(m-3x) -i?oo(0)) e-VdA 

V V ' 

-Am 

+ / /(^) J {|A|>cV^} (^m.( m ~M --Roo(O)) e"VdA. 

•/IS, 

* v ' 

The estimate (2.18a) coupled with the dominated convergence theorem imply that A m — > as m — > 
00. The estimate (2.18b) and the dominated convergence theorem imply that B m — > 00 as m — > 00. 

□ 

Appendix A. Gaussian measures and Gaussian vectors 

For the reader's convenience we survey here a few basic facts about Gaussian measures. For more 
details we refer to [6]. A Gaussian measure on R is a Borel measure 7^, v > 0, m G R, of the form 

lnA x ) = 7y^ e 2v dx - 
The scalar [i is called the mean, while v is called the variance. We allow v to be zero in which case 

Jfj,,o = &n = the Dirac measure on R concentrated at p. 
For a real valued random variable X we write 

XeN(n,v) (A.l) 

if the probability measure of X is 7^. 

Suppose that V is a finite dimensional vector space. A Gaussian measure on V is a Borel measure 
7 on V such that, for any £ 6 V v , the pushforward (7) is a Gaussian measure on R, 

One can show that the map V y 3 £ i-> G R is linear, and thus can be identified with a vector 
V called the barycenter or expectation of 7 that can be alternatively defined by the equality 
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Moreover, there exists a nonnegative definite, symmetric bilinear map 

S:y v xV v ^l such that a(£) 2 = V£ G V v . 

The form S is called the covariance form and can be identified with a linear operator S : V v — > V 
such that 

E(£,t/) = (£,S77>, V^eV v , 

where (—,—) : V v x V — )■ R denotes the natural bilinear pairing between a vector space and its 
dual. The operator S is called the covariance operator and it is explicitly described by the integral 
formula 

(£, Sri) = A(£, rj) = / {£,v - ii~)(t),v - fjt y )dry(v). 
Jv 

The Gaussian measure is said to be nondegenerate if 5] is nondegenerate, and it is called centered if 
/x = 0. A nondegenerate Gaussian measure on V is uniquely determined by its covariance form and 
its barycenter. 

Example A.l. Suppose that U is an n-dimensional Euclidean space with inner product (— , — ). We 
use the inner product to identify U with its dual U y . If A : U — > U is a symmetric, positive definite 
operator, then 

dry A (x) = J- e ~UA-^u, u ) i d - (A 2) 

is a centered Gaussian measure on U with covariance form described by the operator A. □ 

If V is a finite dimensional vector space equipped with a Gaussian measure 7 and L : V — > U is 
a linear map, then the pushforward L*7 is a Gaussian measure on U with barycenter 

and covariance form 

S Lt7 : C/ v x U y -> M, S i , 7 (7 7 , 7 7 ) = S 7 (L V ? ? , J L V ? ? ), Vr? G J7 V , 

where Z v : J7 V — > V v is the dual (transpose) of the linear map Observe that if 7 is nondegenerate 
and Z is surjective, then L*7 is also nondegenerate. 

Suppose (§, /i) is a probability space. A Gaussian random vector on (S, fi) is a (Borel) measurable 
map 

X : § — > V, V finite dimensional vector space 

such that X*^ is a Gaussian measure on V. We will refer to this measure as the associated Gauss- 
ian measure, we denote it by 7x and we denote by Sx (respectively S(X)) its covariance form 
(respectively operator), 

V x (£i,t2) = E((S 1 ,X- E{X) ) (6 , X - E{X) ) ) . 

Note that the expectation of 7x is precisely the expectation of X. The random vector is called 
nondegenerate, respectively centered, if the Gaussian measure 7x is such. 

Let us point out that if X : § — > U is a Gaussian random vector and L : U — > V is a linear map, 
then the random vector LX : S — > V is also Gaussian. Moreover 

E(LX) = LE(X), V LX (Z,0 = Vx(L x/ Z,L x/ 0, V£ € V y , 
where £ v : F v -> t/ v is the linear map dual to L. Equivalently, S'(LX) = LS(X)L V . 



16 



LIVIU I. NICOLAESCU 



Suppose that Xj : S — > V±, j = 1,2, are two centered Gaussian random vectors such that the 
direct sum X\ © X 2 : S — > V\ © V 2 is also a centered Gaussian random vector with associated 
Gaussian measure 

1x^X2 = Px 1 ®x 2 (xi, x 2 )\dxidx 2 \. 

We obtain a bilinear form 

cov(X 1 ,X 2 ) : V\ x V^R, cov(X u X 2 )(£i, 6) = £(6,6), 

called the covariance form. The random vectors X\ and X 2 are independent if and only if they are 
uncorrelated, i.e., 

cov(Xx,X 2 ) = 0. 

We can then identify cov (X\, X 2 ) with a linear operator Cov (X\, X 2 ) : V 2 — > V\, via the equality 

E{£ 1 ,X 1 )(t- 2 ,X 2 ))=cov{X u X 2 ){Z 1 ,Z 2 ) 

= (£ 1 ,Cov(x 1 ,x 2 )tl), V6evy, 6e^> 

where £jj G denotes the vector metric dual to £2- The operator Cov(X\, X 2 ) is called the 
covariance operator of X±, X 2 . 

The conditional random variable (X\\X 2 = x 2 ) has probability density 

f \ _ p Xl ®x 2 (xi,x 2 ) 
P^x^xx) - J ViPXieX2(xuX2) \ dxi \- 

For a measurable function / : V\ — > R the conditional expectation E(f(Xi)\X 2 = x 2 ) is the 
(deterministic) scalar 

E(f(X 1 )\X 2 = x 2 )= / f(xi)p {Xl \ X2=X2) {x l )\dx 1 \. 

JVi 

If X 2 is nondegenerate, the regression formula, [4], implies that the random vector (X1IX2 = x 2 ) is 
a Gaussian vector with covariance operator 

S(X 1 \X 2 =x 2 ) = - Co^^Sp^^Co^Xs,^), (A.3) 

and expectation 

E(X l \x 2 = x 2 ) = Cx 2 , (A.4) 

where C is given by 

C = Cov(X li X 2 )S(X 2 )~ 1 . (A.5) 

Appendix B. A class of random symmetric matrices 

We denote by S m the space of real symmetric m x m matrices. This is an Euclidean space with 
respect to the inner product 

(A, B) := ti(AB). 

This inner product is invariant with respect to the action of SO(rrc) on S m . We set 



Eij, i = j 

i<j. 



Ejn : 



The collection {Eij)i<j is a basis of S m orthonormal with respect to the above inner product. We set 

- J a *i' i = 3 

V2a i? -, i < j. 
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The collection (aij)i<j the orthonormal basis of dual to (Eij). The volume density induced by 
this metric is 

\dA\ := Y[daij = 2s(?) JJ day. 

Throughout the paper we encountered a 2-parameter family of Gaussian probability measures on S m . 
More precisely for any real numbers u, v such that 

v > 0, mil + 2v > 0, 

we denote by 2>m V the space S m equipped with the centered Gaussian measure dT U:V (A) uniquely 
determined by the covariance equalities 

E(aijdke) = u8ij5ki + v(5ikSj£ + SuSjk), VI < .k,£ < m. 

In particular we have 

E(a^) = u + 2v, E(auajj) = u, E(a 2 ij ) = v, VI < i / j < m, 

while all other covariances are trivial. The ensemble S ? ° 7 f is a rescaled version of of the Gaussian 
Orthogonal Ensemble (GOE) and we will refer to it as GOE^. 

For u > the ensemble S^j" can be given an alternate description. More precisely a random 
A G §m V can be described as a sum 

A = B + Xt m , B G GOE^, X G iV(0, u), B and X independent. 

We write this 

§r = GOE v m +N(0,u)l m , (B.l) 
where + indicates a sum of independent variables. 

The Gaussian measure dT u v coincides with the Gaussian measure dT u+ 2 V:UjV defined in [ 5, App. 
B]. We recall a few facts from [ 1 6, App. B]. 

The probability density dT UjV has the explicit description 

dTu,M) = m{m+l ] - e-^-i^ 2 \dA\, 

where 

D(u, v) = (2v) (m " 1)+ (™) (mu + 2v), 

and 

,1/1 1\ u 



m \mu + 2v 2v J 2v(mu + 2v) 
In the special case GOE^ we have u = v! = and 

dT 0jV (A) = Lp— . e -^^ a |dA|. (B.2) 

(2irv) * 

We have a Weyl integration formula [2] which states that if / : § m — > R is a measurable func- 
tion which is invariant under conjugation, then the the value f(A) at A G S m depends only on the 
eigenvalues X\(A) < • • • < X n {A) of A and we have 

if / \ m \ 2 

E GO Ev m ( f(X) ) = — -- / /(A 1; ...,Xm)[ \Xi - Xj\ TT \dXt ■ ■ ■ dX m \, 

Zm( - V > J ^ m \i<r<j<m j t=l 

V v ' 

= '-Qm,v (A) 

(B.3) 
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where the normalization constant Z m (v) is denned by 



i(v) = / \Xi - Xj\T\e~^\d\i- ■ ■ d\ m \ 

jRm i<i<j<™ i=i 

r m x 2 

/ \X i -\ j \T[e-*\d\ 1 ---dX r , 

K Ki<j<m t=l 



The precise value of Z m can be computed via Selberg integrals, [2, Eq. (2.5.11)], and we have 

m t-,/9 \ m 



Z m = (2vr)f ml J] = 2^m\ ]JT £ . (B.4) 
i=i i( 2^ j= i W 

For any positive integer n we define the normalized 1-point corelation function p n>v (x) of GOE^ to 
be 

Pn,v{x) = r . -, r / Qn,t;(») <^2, • • • , A n )(2Ai ■ ■ ■ C?A n . 

For any Borel measurable function / : R — > R we have [7, §4.4] 

-S G OE«(tr/(X)) = / /(A)^(A)dA. (B.5) 
The equality (B.5) characterizes p n) „. Let us observe that for any constant c > 0, if 

A e GOE^ <==>cA G GOE^ . 
Hence for any Borel set B C R we have 



Pn,c 2 v(x)dx= / p n ,v(y)dy. 
cB JB 

We conclude that 

CPn,c 2 v{cy) = Pn,v(y)i Vn > c > ?A ( B - 6 ) 

The behavior of the 1-point correlation function p n>v {x) for n large is described by Wigner semicircle 
law which states that for any v > the sequence of measure on R 

Pn,vn~i{ x )dx = nzp n ^ v {nzx)dx 
converges weakly as n — > oo to the semicircle distribution 

poc,v(x)\dx\ = -f{|x|<2 v M^~V / 4u - x 2 \dx\. 

The expected value of the absolute value of the determinant of of a random A G GOEJj, can be 
expressed neatly in terms of the correlation function Pm+i,v More precisely, we have the following 
result first observed by Y.V. Fyodorov [9] in a context related to ours. 

Lemma B.l. Suppose v > 0. Then for any celwe have 

#GOE^(|det(>4-cl m )|) = 2§(2u) I2 £ J T (j^j^j e^p m +iA c )- 



COMPLEXITY OF RANDOM SMOOTH FUNCTIONS 



19 



Proof. Using the Weyl integration formula we deduce 

-Egoe 



i r x 

v m (\det(A-cl m )\) = — — / T e -^| c -A i |TT|A i -A i |dAi---dA, 

Z m (v) J Rm f± 

f c 2 m *i 

/ TTe~4t|c - Aj| 1 \ \\i - \j\d\i ■ ■ ■ d\ m 



(m 



Z m {v) i=1 
— ^ — TT / ^m+lM ) A l; • • • > A m) a M ' ' ' aA m 

Z m (v) Z m+1 (v) J Rm 

£ <£ 
e*-»Z m+ i(v) m+ieiv Z m+ i 
= ^r-^ Pm+i,„(c) = u 2 - Pm+l,«(c) 

, i—, s m+1 c 2 f m +1\ , . 3 , , m+i / m + 3 \ c 2 

+ l)v/2(2u) — e^r ( — ^— \ p m+1)V (c) = 22(2v)—T ( — ^— J e^ m+ i,„(c) 



□ 



The above result admits the following generalization, [3, Lemma 3.2.3]. 
Lemma B.2. Let u > 0. Then 

_ , . , . . . . N 3 m + 1 / 77i + 3 \ 1 /" . . (c-x) 2 x 2 

E s *,v ( | det(A - ct m )\ ) = 22 {2v) 2 rl — - — j -j== J p m+ltV (c - x)e ^ 2u dx. 
In particular, ifu = 2kv, k < 1 we have 

, /Ty, _1_ Q\ 1 /• 1 ( x I f 2 -)2 j ( f l + 1 ) c2 

23 e »4|det(A-cl ro )|)=23(2i;)^ fcj 4 " di, 

(A := c — 



3 



2,(2,)T|_j-^=y /m+1 ,„(A)e . 



where 



t 2 - 



1-1 1"*' 



Proo/ Recall the equality (B.l) 

§r = GOE^+iV(0,u)l m . 

We deduce that 

E sT ( | det(A - cl m )\ ) = E(det(B + (X - c)l)| ) 
= -tL= / ^goe^ ( I det(J3 - (c - X)t m )\ | X = x)e-£dx 

If 2 
= -= / E G ow m (\det(B - (c- x)l m )\)e~^kdx 
y 2ttu Jr 

a, , m±i /m + 3\ 1 f , . (e-x) 2 x 2 

= 22(2v) 2 r( — - — \-j==J p m+ i, v (c-x)e 2udx. 

Now observe that if u = 2fcu then 

(c x) 2 X 2 X 2 1 , o _ o, 

V . ; - — = -yr- + -p- (a; 2 - 2cx + c 2 
4v 2u Akv 4v 
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1 



( 



-^x 2 - 2cx - c 2 t\ j + 



1 



(x + t 2 k c) 2 + 



c 




4w 



k 



Avt 2 



4/> 



□ 
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